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Abstract 

In the minimal supersymmetric standard model (MSSM) we calculate the one- 
loop finite-temperature corrections from the squarks-Higgs bosons sector, recon- 
struct the effective two-Higgs-doublet potential and study possibilities of the elec- 
troweak phase transition in the full MSSM (m//±, tg/?, At^h, f^, ^q, mu, mo) pa- 
rameter space. At large values of At^t and fi of around 1 TeV, favored indirectly by 
LEP2 and Tevatron data, the threshold finite-temperature corrections from triangle 
and box diagrams with intermediate third generation squarks are very substantial. 
High sensitivity of the low-temperature evolution to the effective two-doublet and 
the MSSM squark sector parameters is observed, but rather extensive regions of 
the full MSSM parameter space allow the first-order electroweak phase transition 
respecting the phenomenological constraints at zero temperature. 
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1 Introduction 



The absence of antimatter in the Universe (the baryon asymmetry), a small ratio of the observed 
number of barions to the observed number of photons n^/n^ ~ 6 x 10^^° and the absence 
of light {niH ~100 GeV) CP-even Higgs boson signal at LEP2 and Tevatron energies lay a 
specific claims to models of particle physics. The baryon asymmetry and an extremely small 
UB/n^ could be understood on the basis of Sakharov conditions, which are respected at the 
electroweak phase transition, expected to take place at the temperature of the order of 10^ 
GeV [1]. The generation of nonzero vacuum expectation value v of the scalar field breaks 
the electoroweak symmetry SU{2)f x U{1)y to the electromagnetic symmetry U{l)em- It is 
well-known [2] that in the simple isoscalar model with the standard-like Higgs potential U (</?) = 
— |//^</7^-|-|A<^^, describing a thermodynamically equilibrium system of the scalar particles at the 
temperature T, the equation for the vacuum expectation value v{T) has two solutions: v(0) = 
and v'^{T) = /U^/A — demonstrating the second order phase transition at the critical 

temperature Tc = \/X = 2v{0), see Fig.la. The thermal Higgs boson mass ml = —jj.'^+XT'^ /A 
vanishes at the critical temperature Tc thus restoring the spontaneously broken symmetry. 



(a) 



(b) 



Figure 1: Contours on the {v,T) plane for (a) the second order phase transition (fraeture point) and 
(b) the first order phase transition at the critical temperature Tc (dashed vertical line). 

However, in the cosmological evolution the stages with thermodynamically non equilibrium 

plasma and the first order phase transitions (see a typical v{T) contour in Fig. lb) arc very 
important, so such simple picture in combination with the standard model CP- violation by 
means of the CKM mixing matrix turns out to be not sufficient to justify the observed ratio of 
baryon number to entropy. The situation becomes better in the minimal super symmetric model 
(MSSM) where sparticles, extended two-doublet Higgs sector with the two background fields 
and nonstandard sources of CP-violation provide a number of new possibilities. In a number of 
approaches [3] the electroweak phase transition is defined by evolution of the finite temperature 
effective Higgs potential involving the cubic term in the background scalar fields Vi,V2. The 
larger this term is, the stronger pronounced turns out to be the first order phase transition, which 
is essential for consistency with the Higgs boson mass beyond the LEP2 exclusion uih <115 
GeV. Enhancement of the cubic term in the MSSM at the one-loop level is substantial in the 
class of MSSM scenarios with a light right stop [4]. Temperature loop corrections from the stop 
and other additional scalar states could be large and lead to the first order phase transition, 
the intensity of the latter depends on ^ = v{Tc)/Tc, where v{Tc) = Jvi{Tc) + V2{Tc) is the 
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vacuum expectation value at the critical temperature Tc. The electroweak baryogenesis could 
be explained if v(T^/Tc > 1 [5], the case of strong first order phase transition. 

In a number of analyses the MSSM finite-temperature effective potential is taken in the 
representation 

V,ff{v, T) = 1^0(^1, ^2, 0) + V^{m{v), 0) + V^{T) + (1) 

where Vq is the tree-level MSSM two-doublet potential at the SUSY scale, Vi is the (non- 
temperature) one-loop resumed Coleman- Weinberg term, dominated by stop and sbottom con- 
tributions, Vi(T) is the one-loop temperature term and Vrmg is the correction of re-summed 
leading infrared contribution from multi-loop ring (or daisy) diagrams. The MSSM relations 
between the SU{2)l x U{1)y gauge couplings g2 and g\, and the quartic parameters Ai,2,3,4 
of the potential Vo(fi,f2,0) are very restrictive. Only two additional parameters tg/3 = V2/V1 
and mH± (charged scalar mass) determine the zero-temperature two-doublet Higgs sector at 
tree-level. The one-loop radiative corrections, both logarithmic and non-logarithmic generated 
at the threshold Msusy-i can change strongly the tree- level picture. They depend on the pa- 
rameters (A,6, A*, '^Q-i fnu-i 'itT'd) of the scalar quarks-Higgs bosons interaction sector. In most 
cases for the analysis in the representation (1) numerical methods are used to find the critical 
temperature Tc, for example, by solving the equation for the determinant of second derivatives 
of the potential (1) at Vi^2 =0 [6]. Then the two background fields fi,2(T'c) are found at the 
minimum using the minimization conditions (i.e. the absence of linear terms of the effective 
potential representation in the "shifted" fields). The first order phase transition strength is 
dependent on the cubic term ETv^ which appears from the infrared region. 

Numerical high-precision Monte Carlo simulations on the lattice [7] have been developed 
and applied to MSSM in connection with the infrared problem [8] inherent to all analyses based 
on the effective potentials. Infrared divergences appear in the integration over bosonic static 
{cuq —0) Matsubara modes, which in the loop expansion for the three-dimensional momentum 
space correspond to the intermediate massless bosons. The non-perturbative investigations of 
the problem have been performed in the framework of high-temperature dimensional reduction 
[9, 10], when an effective three-dimensional MSSM with the same Green's functions as in the 
four- dimensional MSSM for the light bosons is constructed [11, 12, 13] by integrating out pertur- 
batively the non-static modes. The corrections from squarks and gauge bosons are introduced 
after the reduction to the three-dimensional model. 

In order to cover the temperature range from very low temperatures to the temperatures 
of the order of critical, in our analysis we are using different approach developed in [14, 15] 
for the general (non-temperature) two-Higgs doublet potential with complex-valued parameters 
/x^2) -^5) ^6 and A7, which violates the CP-invariance explicitly. However, in this publication we 
consider a simplified situation of the Higgs potential in the CP conserving limit (the imaginary 
parts of the effective parameters As^ej and 2 taken to be zero) . The full MSSM effective 
potential in the generic $1, $2 basis has the form 

C/e//($l, $2) = - /i?($I$l) - f^l{<^l^2) - /iL(*l*2)- 1)12 (*2*l) + Al($I$l)' + A2($^$2)' (2) 
+A3($I$l)($|$2) + A4($I$2)($|$l) + y ($I$2)(M*2) + ^ ($^$l) ($^$l) + 

+A6(<I'I<I>i)($I$2)+ Ae ($l$i)($^$i) + A7(<I'5$2)($I$2)+ At ($5$2)(<^'5<J'i) 

where the background fields (vev's) are = {0,vi)/V2 and ($2) = {0,V2)/V2. The temper- 
ature corrections from squarks, both logarithmic and non-logarithmic (at the SUSY threshold) 
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are incorporated to Ai, ...Ay. In [14, 15] (see also [16]) a nonlinear transformation for masses 
and mixing angles Aj = Xi{a, (3,mh,mH,mA,TnH±\Q, i = 1,...5 to the Higgs bosons mass 
basis can be found for a general case {h,H and A are the neutral and H~^, H~ are the charged 
Higgs bosons, a is the h-H mixing angle, tg/3 = V2/V1) 
22 2 

C/e//($i, $2) =^ ^{hh)+'^{HH)+'^{AA)+mjj±{H+H-)+h, H, A, interaction terms 

Zi 

(3) 

which allows to work with symbolic expressions for the temperature-dependent Higgs boson 

mass eigenstatcs. 

In section 2 we calculate various one-loop temperature corrections to the potential. Section 
3 contains some examples of the electroweak phase transition for the finite-temperature effective 
potential reconstructed in the full MSSM parameter space. The potential of scalar quarks - Higgs 
bosons interaction and some technical details of evaluation can be found in the Appendix. 



2 Finite temperature corrections of squarks 

In the finite temperature field theory Feynman diagrams with boson propagators, containing 
Matsubara frequencies lo^ — I'KnT [n — 0, ±1, ±2, ...), lead to structures of the form 



°o /■ (ik 

/[mi,m2,...,m6] = T ^ / TTT^ 11 77^^ — 2T' (4) 



Here k is the three-dimensional momentum in a system with the temperature T . In the following 

calculations first we perform integration with respect to k and then take the sum, using the 
reduction to three-dimensional theory in the high-temperature limit for zero frequencies. At 
n 7^ the result is [17, 18] 

7[mi,m2,...,m,] =2r(27rrf-^^ ^^^^ ^-^^Y^ 5(M,6-3/2), (5) 

where 

S(M. 6 - 3/2) = I j^^^^,, - (^)^ (6) 

For 6 > 1 the parameter rr? is a linear function dependent on and the variables {dx) of 
Feynman parametrization, which are the integration variables in (6). At the integer values of 
h the integrand in (3) is a generalized Hurwitz zeta-function [19]. Note that for the leading 
threshold corrections to effective parameters of the two-doublet potential 6 > 2, so the wave- 
function renormalization appears in connection with the divergence at 6 = 2 (which is suppressed 
by vertex factors, see [14]). 

A number of integrals can be easily calculated. We calculate the integral 

dk 11 

(27r)3(k2 + a?)(k2 + ai) ~ 47r(ai + a^) ' 

taking a residue in the spherical coordinate system. Here 0^.2 are the sums of squared frequency 
and squared mass, see (4). Derivatives of Jq with respect to a\ and ai can be used for calculation 
of integrals 

^iL0i,02j y (27r)3(k2 + a2)2(k2 + al) 201^01 87101(01 + 02)2' ^ ^ 



Jo[ai, 02] - / - , (7) 
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J2[ai,a2\ I (27r)3(k2 + a?)2(k2 + ai)2 4aia2 9ai9a2 8710102(01 + 02)3" ^ ' 



and ^ 



, dk 1 1 

JsK, a2, 03J = / (27r)3(k2 + o2)(k2 + oi)(k2 + oi) 47r(oi + 02)(oi + 03)(o2 + 03)' 



,r 1 I 2ai + 02 + 03 . s 

J4Lai,02,03j- j (27r)3(k2 + af)2(k2 + oi)(k2 + ai) ~ 87roi(oi + 02)2(01 + 03)^(02 + 03) ' ^ ^ 



Thus, the procedure of Feynman paramctrization is not used. Substituting af = An'^n'^T'^ + ml 
and O2 = 47r2n2T2 + m| to (7) and taking the sum over Matsubara frequencies after the inte- 
gration we get 

00 CO -j^ 

/o[mi,m2]= V Jo[mi,m2]= V , = , (12) 

n=—cx),n^O n=— oo,n^0 



or, after redefinition of mass parameters M1.2 = mi.2/27rT the temperature corrections to effec- 
tive potential are expressed by summed integrals 

1 °° 1 

7i[Mi,M2l = — y , , (13) 

6471^^' n=-t:n^o ^Ml + n'i^Ml + n2 + y'Mf + n^f 

I 00 I 

/2fMi,M2l = ^ V ^ , (14) 

2567r5r4 yM^:^^ /](^|^( /^i^^:^^ Jmf+^f 



Note that the series (12) are divergent, but the derivatives (13) and (14) are convergent for all 
M1.2. In the following it will be convenient to keep separately the terms for zero and nonzero 
modes in the sum. Both terms will be temperature-dependent since the zero-mode integrals 
coincide with (7)-(9), where of = mf and the factor T should be accounted for. Numerical 
check of the zero temperature limiting case T — > demonstrates that the non-temperature field 
theory results are successfully reproduced. In the high-temperature limit the zero mode gives 
dominant contribution in agreement with a known suppression of quantum effects at increasing 
temperatures. 

The sum of integrals (13) and (14) can be expressed by means of the generalized zeta- 
function. Such forms can be derived if we introduce Feynman parameters in the integrand of 
(7) 

' (15) 



[k2 + m2][k2 + m2] h ([k2 + m2]x + [k2 + m2](l-x))2' 
and redefine k — ^ p = k/27rr, M2(M„, M^, x) = {Ml - Ml)x + M^. Then we get 

1 I dx 



I-kTyL (p2 + n2 + M2)2- ^^^^ 



[k2 + m2][k2 + m2] (27rT)4Jo (p^ + n^ + M^)^ 
and divergent series for (7) (dk = {2-nTYdp) 

1 /-i °° r dxi 1 

/„[M„.Md = —I E / 7^ (p.^„.^^pp . (17) 



n=— oo,n^0 ■ 



^The same results for J3 and J4 can be found in [11] and [12], where they appear in the context of high 
temperature dimensional reduction. 
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With the help of dimensional regularization or differentiating the integral 

J (27r)3 (p2 + M2) 471 ^T^' 
over the parameter M, the equation (17) can be reduced to 

where ({u,s,t) is the generalized Hurwitz zeta-function [19]: ^ 

oo 2 

C(u, s,t) — Y^ — . 



(18) 



(19) 



(20) 



So in the case under consideration the sums of integrals (13) and (14) can be calculated by 
differentiation of (19) with respect to mass parameters participating in M = M{Ma, Mh,x). 
Differentiation increases the power s in the denominator of (19) giving convergent integrals 



1 d 
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The integrals (21) and (22) are equal to the series (13) and (14), respectively. 



(21) 



-/,) = [\xxil-x)C,{2,-,M^). (22) 






Figure 2: Threshold corrections (left and central diagram) and diagram contributing to the wave- 
function renormalization (right). 



Threshold corrections from the triangle and box diagrams, shown in Fig.2, are 
denoted by AA*'*, They contribute additively to the parameters Aj = Af^"^^ — AA*'*. 

In the following we are using the normalization conventions from [14]. Calculation of the finite- 
temperature diagrams for the general case of complex-valued fi and At^h gives the result (see 
details in the Appendix) 

AAf = ?,hl\^x\^h[mQ, mu] + 3/i^|A|%[mQ, mij]+ (23) 

^^-^^h[mQ, mu] + 2glh[mu, tuq]) 
^Note that (non-generalized) Hurwitz zeta-function is defined by C,{s,t) = X^^o \n+t\' " 
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AX't = :^ht\At\%[mQ,mu]+?>htWl2[mQ,mD\+ (24) 



+(l/^r^%^ + |Ar|)/i[m^,mQ])+ 



+(i/^r^%^ + i^.p|Ki[^i.,^Q])+ 

+/i?|/x|Vtl^^2[mQ, mu\ + /i^lAtHAbl^/ahQ, 771^] + 
+h'^hl{2{AtAb - \n\'^)l3[mQ, mu, mo] + (Ia*!^ + l^ri^lbT - '2AtAb\ij\'^)h[mQ, mu, mo] 

AAf = 6hM'\At\^l2[mQ,mu] + 6ht\fi\^\A,\^h[mQ,mD]+ (26) 

+^(1^1'^!' - |/^r(^i' - Qhl))h[mn, tuq]) - AXf 

AAf^ = 3/iyA^72[mQ, mt;] + 3htfJ,^Alh[mQ, m^] (27) 
AA*'^^ = -3htfxAM''l2[mQ,mu] - 3htfiA,\A,\^l2[mQ,mD]+ (28) 

+^tMt(^^^-^^-^i["^Q,m[/] - gpi[mu,mQ])+ 

+%M6( ^ h[mQ,mD\ ^ h[mD,mQ\) 

AXf' = Ar/2[mQ, mc,] - 3htfiA,\fx\''h[mQ, m^] (29) 

+hliiAb{- ^'^ ^^^^ /i[mQ,mj,] - y/i[m£,, mQ])+ 

+htl^At{ — — ^ —h[mQ,mu] - {3h^ - gl)Ii[mu,mQ]) 

where gi, g2 arc U(l) and SU(2) gauge couplings, fi is the Higgs supcrficld mass parameter, 
At, Ah are the trihnear squarks- Higgs bosons parameters, ht,hb are the Yukawa couphngs and 
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Figure 3: "Fish" diagrams 



mQ,mu,mD denote the scalar quark mass parameters, in terms of which the physical masses 
are expressed. 

Corrections of "fish" diagrams, see Fig. 3, give the following contributions to the effective 
parameters 



-AAf = 



2 9t 
6 



gf 



{I(mQ) + I{mn)) + fl{mu), 



I(mQ) + [h',-^-inmu) + ^I{mn), 



-(AA3 + AA4)^ = ^ {-gt + Hhl - hl)gl - - 2{hl + hl)gl)) I{mQ)+ 
-AA{ = ^ {-gt + 6{hl - hM + 9 {gl -'2{hl + h^) gl + ^hlhf)) I{mQ)+ 



72 



-AA{ = (hi - f )(f - h^Mmg) - h'Xiimu, mo). 
The three-dimensional integrals in (30)- (34) are 

1 



J {mi) 



J{mu,mD) 



4:n{mu + mo) 



(30) 
(31) 



(32) 



(33) 
(34) 

(35) 



see (7), leading to series analogously to (12) and (17). 

The logarithmic corrections for non-degenerate squark masses can be defined following 
[20] and [21]. Schematically in the results of [15]) we replace In [^^^) by In (^^^) : 



AAl''^ = -3g^ (ll^^ - 2>^hlgl + 9 (gt - Ahlgl + l^hf)) In 



AA^^ = (44^1 - 144/.?^? + 36^2^ + 576/i^ - UAglhl) In 



( mQmu\ 



/ mQmu\ 



mi 



(36) 
(37) 



3847r2 



(-11^^ + 18 (hl + h^) gl+ 



(38) 



+9 {gl -2{hl + h^) gl + l^hlhf)) In , 



647r2 



{gt - 2 {hi + hi) gl 



Shlhl 



In 



/ mQmu\ 



(39) 



Large logarithms not connected with the renormahzation group appear also in the wave-function 
renormalization yield, see below. 

It is known that in order to renormalize the the X^p'^ theory, one needs to renormalize the 
self-coupling and the mass of the scalar field. If the A(/7^ theory is supplemented by fermions 
with interactions defined by the Yukawa term, an additional wave-function renormalization 
is necessary. Similar situation takes place in the two-doublet model. Expanding the self- 
energy diagram (see the insertion to the leg in Fig.2, right) calculated with non- degenerate 
masses at finite temperature, we get at =0 the wave-function renormalization (w.f.r.) 
correction, which is defined by a factor in front of p^. At zero temperature two ways of w.f.r. 
calculation can be used [23]. Our calculation is based on the integration of convergent w.f.r. 
contribution over the momentum squared, previously which has been used in differentiation. The 
standard subtraction scheme at zero momentum (BPSZ-scheme) in the divergent expression 
for the self-energy contribution, when the divergent pole part is subtracted, turns out to be 
not convenient at finite temperatures, because in summation over Matsubara frequencies not 
divergent integrals, but divergent series must be subtracted. Following [14] we can write 



A Xt -\i9l + 9l)Ki. A Xt = \{gl + gl)A!,,, 



(40) 



^Xt^-\{9l-9l){A,+^,,), i^Xt^-\9l{Ki + ^22), AA^^^^O, 
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{gl + 9l){K2 - A',:) = 0, A A^*-- = -igl + gl){A',, - A',, 



0. 



where A matrices 2x2 are ^ 



_ 2^ 

^'^-^ 24 TT 



F{mQ,ml,T) 



-jiAu 



-li*A}j 



+ ^D,A^^x)}{l--l), (41) 



include the series (compare with Eq.(113) in [9], taking into account differentiation to get the 
finite w.f.r. yield) 



+00 



F{mlmlT) = T 



1 



T 



-oo {^ml + (27rnr)2 + ^ml + {2T:nTff 

1 



(42) 



+ 0O 



{nil + 1^2)'^ 



1 {J ml + (27rnT)2 + J ml + (27rnT) 



\2\3 



The sum of all w.f.r. corrections to X^^j vanishes. 

It is useful to check that the finite temperature corrections are reduced to the structures of 
zero-temperature MSSM, which play a role of boundary condition at T—0. Indeed in the 



^The equations in [14] are given for the general case of complex-valued /z, A^^h- 
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limiting case of T—0 and degenerate squark mass parameters all equal to Msusy the threshold 
corrections given by Eq.(23)-(29) are reduced to previous zero-temperature results [14, 22]. For 
example, let us take AAi for = nib = Msusy 

AAi = Shtlfil^'UMsusv] + ^ht\A\^l2[MsusY]+ (43) 

+hl\A\'C-^^^t:f^h[MsusY] + m - gDhWsusY]), 
where the integrals are 

r d'^k 1 11 

h\MsvSY] = -J p^j, 1^1^, ^ ^,^2^^^_j3 = -16^ mq^' 

hlMsusr] ^ I ^(P + m|,«k)* = TS^si; 
Transformation to Minkowski space leads to the change of sign in (44). The equality of the 
temperature series for /i 2 to the symbolic expressions for the integrals can be numerically 
verified. 

In the limiting case of T=0 and different squark mass parameters the reduction of (13) and 
(14) to the four-dimensional Ii and I2 can be achieved using 

1 1 d 1 

= , (46) 

[p^ -|- m^]^[p^ -|- ml] 2ma drUa + ?7T.^]b^ + i^t] ' 

Differentiating (46) with respect to 

(47) 



jp2 _|_ j^2j2jp2 _|_ j^2j2 duib [p^ -l- m^]^[p^ -|- mf\ 

then using Feynman paramctrization (15). differentiation in the same way as in (46) and (47), 
and dimensional regularization to integrate over the four-momentum p with the following inte- 
gration over the Feynman parameter, we arrive at 



' ' (27r)4[p2 + ^2]2[p2 + ^2] lQT,\ml - ml^^ ' ^ ' 



d'^p 1 mi - ml - {ml + ml)ln— 

(27r)^ [p2 m^]2[p2 _^ ^2j2 87r2(m^ — m^)3 



rrn, 



(49) 



In the limit ma = mi, these formulas coincide with the expressions for degenerate squark 
masses (44) and (45). 

In calculations of the temperature dependent parameters Xi{T) of the effective MSSM po- 
tential at moderate temperatures we used truncated series with fifty terms (50 Matsubara 
frequencies). Relative contributions of the remaining terms are less than 10~^ percent at T=50 
GeV, decreasing with an increasing T. At small temperatm^es of the order of a few GeV an 
acceptable accuracy is achieved with 1000 terms. The effective parameters Xi{T) are less than 
one, justifying the perturbative approach, as a rule, at the squark mass parameters around 
several hundred GeV. However, strong parametric dependence is observed here, for example, at 
the squark mass parameters 200, 500 and 800 GeV the criteria Xi{T) <1 is valid up to T ~ 860 
GeV, while taking degenerate squark masses at 600 GeV we found that at T > 600 GeV the 
perturbative regime cannot be used. 
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Figure 4: The zero-temperature two-doublet Higgs potential Uo{vi,V2), see (2), at the scale MsusY- 

3 Thermal evolution and the critical temperature 

In view of the effective two-doublet potential structure defined by (2) one could assume that 
broken symmetry of Ueff{vi,V2) with a local minima at T =0, Vi^2 7^0 appears in the sum of 
the potential terms with fif, /ig and /ifg of dimension 2 in the fields, which form a 'saddle' 
(a hyperboloid in the (^1,^2) space), and of the dimension 4 terms Ai^...7 which are increasing 
quartically, being unbounded from above. However, the situation is more involved because fil 2 12 
and Aj respect a number of constraints. In this section we are going to describe roughly some 
possible scenarios of temperature evolution in the effective two-doublet MSSM Higgs sector with 
threshold, logarithmic and wave-function renormalization one-loop corrections. The squark mass 
parameters in the following numerical calculations are rriQ =500 GeV, =200 GeV, mn =800 
GeV, giving substantial mass splitting. With these parameters the third generation squark 
eigenstates rritj^^ and mt^ ^ have masses in the range from 400 GeV (mbj to 900 GeV (mt^) at 
the values of At^ and /i of the order of 1 TeV (only an extremely high-tangent (3 ~50 region is 
sometimes not suitable). Large difference of the stop masses is necessary to respect constraints 
following from the LEP2 experimental limit ruh >115 GeV. The values of tg(3 above 5 and 
large soft supersymmetry breaking parameters At^f, and ^ of the order of tuq also lead to an 
acceptable Higgs boson mass rrih (but weaken the strength of the electroweak phase transition 
if taken too large). At the same time substantial threshold corrections appear in the MSSM 
scenarios with large At^b and /i, like the BGX scenario [24] and the CPX scenario [25], or the 
regions of MSSM parameter space close to BGX and CPX. 

The potential (2) includes temperature-dependent parameters Aj(T), i=l,...7, and Vi^2iT), 
see Eq.(23)-(29), which define the thermal evolution from some high temperature T of the 
order of several hundred GeV down to zero. At the SUSY scale, where Ai = A2 = {gf + 92) /8, 
A3 = {g2— gf)/4:, A4 = —(72/2 and A5 = Ae = A7 = 0, the zero-temperature potential Uo{vi,V2) = 
— {gf + 5'2)('^i ~ "^2)^/32 is unbounded from below and has two 'fiat directions' Vi = ±V2, see 
Fig. 4. Threshold corrections at zero temperature can be found in [22]. As a rule they transform 
the decreasing function in Fig.4 to a saddle configuration, slowly increasing along one of the 
'flat directions' and more rapidly decreasing along the other. Necessary ingredients of the two- 
doublet Higgs system are the minimization conditions which set to zero the linear terms in the 
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-200 



-200 



Figure 5: Development of the saddle configuration in the potential U(.ff{vi,V2), see (2), at the critical 
temperature Tc =120 GeV. The squark sector parameter values are mq =500 GeV, mij =200 GeV, 
mo =800 GeV, At = Ab =1200 GeV, /i =500 GeV, the charged Higgs boson mass mH± =150 GeV. 
Horizontal plane corresponds to t^e// =0- 

physical fields h, H and A and ensure a local minimum at any point of the surface Ueff{vi,V2) 
in the background fields space (see e.g. [14]) ^ 

fil = Ait;^ + (A3 + A4 + ReA5)y-Re/i?2tg/3+^(3ReA6Ctg/3 + ReA7tg/3), (50) 
= X2vi + (A3 + A4 + ReA5)y - Refil^ctgp + ^(ReAectg/? + SReArtg/?), (51) 
where (m^i is the CP-odd scalar mass) 

Re/ii2 = sin/?cos/3[m^ + y(2ReA5 + ReAectg/? + ReAytg/?)]. 

Important input parameters of the two-doublet potential are tg/5 = V2/V1 and the charged Higgs 
boson mass 

"^H± ="^w^ + "^A-y(ReAA5- AA4) (52) 

where the effective temperature-dependent mass of the longitudinal ly-boson is m^^(f,T) = 
m^(t') + nvi/^(T), IIwl{T) = bg^T"^ /2 (with the one-loop Standard Model and third-generation 
squarks contributions included in the polarization operator; = v'^gl/2). We denote A345 = 
A3 + A4 + ReAs. If in the process of thermal evolution, when the system moves along some 
trajectory in the vi{T),V2{T) plane, we require the minimization of U with respect to the scalar 
fields oscillation in the extremum vi{T)^V2{T) and continuously admit the interpretation of the 
system in terms of scalar states h,H and A, then ix\, ix\ and /i^g can be expressed by means 
of the effective parameters Ai ...7 [14]. ^ Only ii\ and /x| are dependent on the direction in the 
(wi,f2) plane, while /i^2 ^^^d Ai,...7 are not. 

Although only the CP-conserving limit is considered, we keep the notation of real parts for the variables 
where a phase factor could appear in the general case. 

^The normalization of A1.2 in [14] is different from [15] by a factor of 2 
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A, GeV A, GeV 

Figure 6: Contours of negatively defined Ai (left, dark grey area) and A2 (right, dark grey area) in the 
{Af = Af,, jj) plane at the temperature 150 GeV. The squark sector parameter values mq =500 GeV, 
mu =200 GeV, mo =800 GeV. 




200 400 600 800 1000 1200 1400 



A, GeV 

Figure 7: Contour of negatively defined determinant A1A2 — A345/4 (dark grey area) in the {A^ = A^, 
/i) plane at the temperature 150 GeV. The squark sector parameter values mq =500 GeV, mu =200 
GeV, mo =800 GeV. 

First we consider the simplified case Ae = A7 = 0. The two-doublet Higgs potential without 
\q and A7 terms has been considered in the context of discrete Peccei-Quinn symmetry [26]. 
In this case the effective potential Ueff{vi,V2) = — (Ai^;^ + \2V2 + A345vfv|)/4 is positively 
defined and unbounded from above if the Sylvester's criteria for the quadratic form C4//(vi, vl) 
is respected 

Ai<0, A2<0, AiA2-^<0 (53) 

At the critical temperature defined by the equation A1A2 — AI45/4 =0 the positively defined 
potential starts to develop the saddle configuration which is unbounded from below, see Fig.5. 
The regions of positively and negatively defined Ai and A2 and the contour for Sylvester's criteria 
(53) are shown in Figs. 6 and 7 at the temperature T =150 GeV in the {A = = Af,, fi) 
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plane. The squark mass parameters mg, mu and rrif) are fixed as mentioned in the beginning 
of the section, the parameters are chosen in the vicinity of the contours which separate 

positively and negatively defined A-parameters in (53). The critical temperature in this case is 
slightly above 120 GeV, insignificantly dependent on the values of {At^, n) if they are changing 
along the right segment of the contour in Fig. 7, separating the light grey and the dark grey 
areas. For the general case of nonzero Ag and Ay defined by Eqs.(28) and (29) the effective 
potential Ueff{vi,V2) = — (Aif ^ + A2f| + A345f^f| + 2A6f?f2 + 2A7f if|)/4 always demonstrates a 
saddle configuration which slopes become steeper with an increase of the temperature. Typical 
shape of Ueff{vi,V2) is shown in Fig. 8. 




Figure 8: The potential Ueff{vi,V2), see (2), with nonzero Xq and Ay at the temperature T =120 GeV. 
The squark sector parameter values mq =500 GeV, mu =200 GeV, mo =800 GeV, At = Ai, =1500 
GeV, fj, =1000 GeV, charged Higgs boson mass =150 GeV. Horizontal plane corresponds to 

Ueff =0. 

So far we assumed that we are working in the framework of effective field theory at the rritop 
energy scale, when the contributions from squarks decouple or a contribution of the potential 
terms with squarks (see the Appendix) is practically constant. However, if it is not the case and 
the Higgs bosons-squarks quartic term is positive definite with the global minimum at the origin 
Vi = V2 =0, the phase transition may occur due to the development of a saddle configuration by 
the fii,fi2 and /if 2 terms. Such situation may take place when the vacuum expectation values 
of charged and colored superpartners participate in the full MSSM scalar potential, possibly 
giving charge and color breaking minima [27]. It is convenient to use the polar coordinates 
fi(T) = v{T) cos (3(T), V2(T) = v(T) sm/3(T) for the vacuum expectation values. The mass 
term of the two-doublet potential has the form 

— v'^ — — — 

Umassiy.P) = - — {fil cos^ P + ulsin^ P) - y/i?2sin2/3 (54) 

By definition at the critical temperature the gradient of Umass{v, (3) is zero along some direction 
in the (fi,f2) plane, then dUmass/dv = and 1/v dUmass/dP = 0; it follows from these two 
equations 

tg2p = {f^V2 - f^'u) M + ^f^'u] = (55) 
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The phase transition is characterized by the critical angle (3{T) which defines the fiat direction ^ 
for the mass term at the temperature in the background fields plane {vi,V2), and at the real- 
valued yUi, /i2 and /ii2 the critical temperature is defined by the equation 



/^?/^2 = /^12 (56) 

For a fixed set of the squark sector parameters ^ the thermodynamical evolution of the effective 
potential is described by a gradUeff =0 trajectory in the three-dimensional (f , T, tg/5) space, 
which is defined by the intersection of the two surfaces, corresponding to the equation (55) for 
the critical angle ("/9-surface"), and the equation (56) for the /ii, fi2 and fiu ("/i-surface"). We 
show the cross sections of /i-surface (calculated without any approximations numerically) by 
the plane at fixed T=0, giving the (f , tg/?) contour, and the cross section at tgP =1 giving the 
(f , T) contour in Fig. 9. The presence of nonzero effective parameters Xq and Ay is essential to 




50 100 150 200 50 100 150 200 250 300 



T,GeV V, GeV 

Figure 9: Cross sections of the /x-surface, see (55), at the temperature close to zero (right panel) 
and tg/3 =1 (left panel, see also Fig.l). White area in the right plot corresponds to the parameter 
values when the effective mass term (54) has a saddle configuration, Af^i, =1800 GeV, fi =2000 GeV. 
The /3-surface is very close to the //-surface. The squark sector parameter values are mq =500 GeV, 
mu =200 GeV, rriD =800 GeV, the charged Higgs boson mass mfj± =150 GeV. 



get the critical temperature of the order of 100 GeV. Useful analytical approximation can be 
obtained using the minimization conditions (50)-(51), then the critical angle /3(T) defined by 
(55) can be expressed as ^ 

tg2P = tg2(3—^ . (57) 

Q, ) 2Ai cos-'/3-2A2 \ i ^"^A i 

^2^ "1^ ^^^2^ A345 + — + 02 



^Flat directions may exist also in the quartic term separately taken, see e.g. [28]. 
Mathematica package [29] with encoded representations of Xi{T) by means of series with n—50 was used to 
scan the MSSM parameter space. At low temperatures the convergence of Matsubara series becomes worse, so 
the number of terms up to n=1000 is needed to reach an acceptable accuracy. 

^In different context analogous relation between tg/3 and tg/3 can be found in [6], where the mass term of 
the form v'^f{(3,T) with f{(},T) = a{(3)T'^ — b{f3) was analyzed for a special case of degenerate squark masses, 
A = ^ = and with the high-temperature expansion. Our quartic potential is very different from the tree-level 
(Si + 92) /8 ^tv'^ cos^ 2f] plus a logarithmic term [6], so the approximate expression 3E/Xt >1 of Bochkarev- 
Shaposhnikov criteria v{Tc)/Tc >1 for the absence of sphaleron in the broken phase is not suitable with nonzero 
threshold corrections. 
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where 

tti = Y + \{Kctg(3 + Avtg/?), a2 = K{tg2(3 - ctg/?) - M{tg(5 + tg2/5). (58) 

If wc assume that /9(T) = /3(0), i.e. only the modulo of v{T) but not the direction in (■^1,^2) 
plane are changed in the process of thermal evolution, then for the critical angle (57) we get 

-^(2A5 + Ae Ctg/3 + A, tg/3) + ^ [ 2A.-2A.tg^/3+Ae(3tg_^-^y)+A.(tg3;3-3tg/3) _ ^^^^ ^ Q_ (59) 

This approximation may be too rough at small mH± , as pointed out in [6] . The saddle configu- 
ration changes not only the shape, but also the horizontal orientation in the process of thermal 
evolution. For the case A5 = Ae = Ay = (59) is reduced to 

^s'-^-lr^ (60) 

ZA2 — A345 

Combining (51) and (59), where we keep only the leading power terms in Aj and omit A5 << 
m^/v^, the equation (56) can be written in the form 

Ai (2A2 - A345)' + A2 (2Ai - A345)' + A345(2Ai - A345)(2A2 - A345) = (61) 

The vacuum expectation value v and mass tua do not explicitly participate in this equation, 
only the dimensionless effective parameters Aj of the quartic potential terms. The left-hand 
side of (61) approaches zero from below as v increases, demonstrating however no solution for 
the saddle configuration. This can be understood quahtatively if we rewrite (61) in the form 
Aictg^/3-|- A2tg^/3-|-2 A345 =0 where the numerical values in the A-pattern, see Fig. 16, calculated 
in the BGX scenario Ai < 0, A2 > and A345 < 0, so in (60) tg/3 < 1. 

Turning back to the case of effective field theory when the squarks decouple at the rritop 
energy scale, the evaluation of thermal masses of Higgs bosons, mixing angles and couplings can 
be done using results of [14]. For example, the thermal evolution of the CP-even Higgs bosons 
h and H is expressed by 

T^l = 4_^m^ + i;2(2Ai4c2 +2A24s2 _2(A3-FA4)c„c^s«s^ + ReA5(44 + 4c2) (62) 

-2c„+/3(ReA6SaC/3 - ReAyCaS/j)), 

m]j = 4-/3"^A + ^^(2Ai44 + 2A244 + 2(-^3 + A4)caC;3SQ,S/3 + ReA5(4sJ + 44) (^3) 

-F2sQ,+^(ReA6CaC/3 ReAySaS/j)), 

where the mixing angle a of the CP-even states h and H is 

_ S2pm\ - t;^((A3 + \a)s2p + 24ReA6 + 2s^ReA7) 

C2i3'm\ - ^^(2Ai4 ~ 2A2S^ - ReA5C2/3 + (ReAe - ?^e\j)s2p) ' 

We show the regions of the (f ,T) plane where the CP-even Higgs boson masses nih and 
iriH are positively defined in Fig. 10 - 12 (shaded areas). "Tachyonic" areas (shown in white 
colour) correspond to the negative squared masses of ruh or mH-i see (3), when the fiuctuations 
of physical fields h and H near the unstable local extremum (wi(T), ^2(7")) grow exponentially 
with time. Configuration of the effective potential Ueff{h,H,A,T) expressed in the physical 
fields h, H in these areas is a saddle or a function with negative or indefinite sign values 
unbounded from below. The 'saddle' temperature Tc —120 GeV (see Fig. 5) is close to the 
temperature, when the thermal mass mh{T) vanishes, only at the low tg/3 values. The heavy 
scalar mass muiT) vanishes at the temperatures substantially higher than Tc- In the scenario 
under consideration at higher tgf3 one should increase the charged scalar mass to respect the 
zero-temperature condition v{0) =246 GeV. 
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Figure 10: In the shaded areas m/i (left panel) and rriH (right panel) are positively defined at the 
parameter values tg/3 =5, mjj± =180 GeV, Atfi =1200 GeV, fi =500 GeV. Isocontours of constant nih 
and TJiH masses are indicated. The squark sector parameter values are uiq =500 GeV, mu =200 GeV, 
ruD =800 GeV. 




50 100 150 200 50 100 150 200 250 300 



T,GeY T,GeY 
Figure 11: The same contours as in Fig.lO at tg/3 =15, mH± =230 GeV 

4 Summary 

Our analysis of the effective MSSM finite-temperature potential is based on a calculation of 
various one-loop temperature corrections from the squark-Higgs sector for the case of nonzero 
trilinear parameters At, and Higgs superfield parameter /i. Quantum corrections are incorpo- 
rated in the parameters Ai^...7 of the effective two-doublet potential (2), which is then explicitly 
rewritten in terms of the Higgs boson mass eigenstates, using the approach developed in [14, 15]. 
The effective parameters Ai(T),...A7(T) include the threshold corrections from triangle, box and 
"fish" diagrams together with the logarithmic and the wave- function renormalization terms. 
The dominant contribution comes from the triangle and box graphs (Fig.2, left and central) 
and can be written in a compact form by means of the generalized Hurwitz zeta-function. 
Temperature evolution of the potential Ueff{vi,V2) expressed in terms of the background 
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Figure 12: The same contours as in Fig. 10 at tg/3 =40, mH± =260 GeV 

fields (fi,f2) is very sensitive to the MSSM scenario under consideration. We are using the 
scenarios with large At^b and /i (about/of the order of 1 TeV), favored by available experimental 
data. When the non-temperature threshold corrections [22] to the tree-level zero temperature 
MSSM Higgs potential, see Fig. 4, are included at some values of mg, m[/, m^, At^ and /i, it 
changes the shape, becoming sign indefinite and unbounded from above, developing a saddle 
configuration. The evolution of a saddle with the increase of temperature is strongly dependent 
on the availability of Ag and A7 terms and the magnitude of possible contributions from the 
MSSM squarks-Higgs bosons sector. If the latter are small (or practically constant in the 
background field space), heavy superpartners of quarks decouple at the top mass scale and 
Ueff{vi,V2) is adequate to the temperature range of the order of 10^ GeV, then the phase 
transition at some critical temperature of the order of 100 GeV takes place in an extensive 
regions of the MSSM parameter space and is described by a development of the saddle structure 
which is defined by the quartic dimensionless parameters Aj(T), i = 1,...5, see Fig. 5. The 
potential Ueff{h,H,A) written in terms of physical Higgs fields (i.e. Higgs mass eigenstates 
h,H and A) demonstrates the spectrum of scalars with positively defined masses which are 
reaching zero at different temperatures, see Figs. 10 - 12. These temperatures rule, 
not close to the 'critical' temperature Tc, when the potential Ueff{vi,V2), Xq = A7 =0, forms 
a horizontal 'narrow gully', so not only the first, but also the second derivatives are zero in 
some direction. It is interesting that the transformation of the potential configuration at zero 
temperature unbounded from below, see Fig. 4, to the saddle configuration unbounded from 
above, see 5, in the framework of large f, and /i MSSM scenarios is rapid if thermal masses 
of scalars do not change substantially. The isocontours for CP-even scalar masses rrih and ttlh 
fall down in an extremely narrow temperature region near T =0, see Figs. 10 - 12, so during 
the 'overturn' of the potential a nearly step-like decrease of v{T) must happen to keep constant 
masses. If the contribution from the squarks-Higgs bosons sector is large and the decoupling 
approximation for squarks is not valid, the phase transition is generated by a saddle configuration 
of the mass term with fil, ii\ and /i^g; dimension 2 with respect to the fields. Although our 
evaluation is performed at the one-loop only and not all possible corrections are considered, 
usually it is possible to adjust yU, At^b (or Xt^b), ^q, rnu and rriD of the MSSM parameter space 
in such a way that the boundary condition for zero temperature f (0) =246 GeV is respected, the 
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critical temperature is of the order of 100 GeV or higher, and the phase transition is of the first 
order. Independently on the temperature evolution scenarios, one should not forget that the 
value of tg/9 at zero temperature must be consistent with the range from 5 to 50-60, provided 
by phenomenological restrictions from LEP2 and Tevatron data for the reactions e~^e~ — > hZ, 
h bb, and pp — > , i — > H^b. In the nearest future useful information about the allowed 
regions of the MSSM parameter space could be provided by the LHC Higgs physics program 
[30]. Availability of the criteria v{Tc)/Tc ~1 for the absence of sphaleron in the broken phase 
deserves a more careful study with the evaluation of radiative corrections from other sources, 
especially the infrared ones. 

We considered only the case of real- valued MSSM parameters At^b and fj,. Generalization 
to the complex-valued parameters (the case of explicit CP violation in the squark-Higgs and 
the two-doublet Higgs sectors) is straightforward with radiation corrections defined by Eqs. (27)- 
(29), where phases of At^b and /i can be introduced. Complex X^^qj lead to the mixing of CP-even 
h,H and CP-odd A scalars resulting in the Higgs bosons without a definite CP-parity /ii,/i2 and 
hs with specific properties, modifying the qualitative picture described above. 

Acknowledgements 

M.D.(MSU) is grateful to Mikhail Shaposhnikov for useful discussion. The work of M.D.(MSU) 
was partially supported by grant NS 1456.2008.2. The work of M.D.(SSU) was partially sup- 
ported by grant ADTP 3341. 



18 



5 Appendix 

Inputs and some details of various quantum corrections calculation, see Eq.(23) and the following 
formulas, are given below. Supersymmetric potential of the Higgs bosons - third generation of 
scalar quarks interaction has the form [20] 



where 



V0 = Vm + Vf + Va + Vq, (65) 

Vm = -fJ-lM^j + M2 (g tg j + M§U*U + M^D*D , (66) 

Vr = rf ($tg) D + Tf aag) U+ Tf (Q t$,) D*- Tf (tQ ^a^^*) W , (67) 
Va = Ag ($1$,) ($1$,) + ($1$,) [Ag (g^Q) + A^jU*U + ASD*d] + (68) 

+Ag (*lg) (g + 1 [Ae,,- (72*,) ^*t7 + h.c] , i,j,k, I = 1,2, 

Vg denotes the four scalar quarks interaction terms, (T2 = ^ ^ ^ , and A^ are defined by 
the tree-level equalities 

A« = diag{i(^2' - glYQ), hi - \{gl - glYQ)}, 

7^ = d±^g{hl - \gl, \gl - hi}, A^ = diB.g{-\glYu, hi + \giY^}, 

A^ = diag{/i^ - \g\YD. A = -hyhj,. 

Here the squark hypercharges are Yq^ = l/3(— 1), Y^^ = 2/3(2), Yu. = —4/3, and the Yukawa 
couplings for the third generation squarks ht = -^^f, hb = • In the general case of 

complex- valued parameters 

r\i.2} ^hu{-fi*;Au}, r(i.^y = hD{AD;-fi*}, (69) 

In order to calculate, e.g., the one-loop threshold corrections (23)-(29) first we extract from the 
potential (65) the triple and quartic interactions presented in Fig. 13 and Fig. 14. The triangle 
and box diagrams which contribute, for example, to the threshold corrections included in Ai 
are shown in Fig. 15 together with symbolic expressions for the temperature one-loop integrals. 
Their sum multiplied by the color factor 3 gives A^'"", see (23). 




Figure 13: The triple interactions extracted from the squark-Higgs sector, see (65). 




Figure 14: The quartic interactions extracted from the squark-Higgs sector, see (65). 
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Figure 16: Histograms for temperature-dependent Aj {i=l,...7) with various quantum corrections in 
the framework of the CPX-hke scenario [25], At = =1000 GeV, ^ =2000 GeV, in the cases of (a) 
degenerate squark masses mg = mt = nib = ^SUSY =500 GeV, zero temperature, (b) degenerate 
squark masses uiq = mt = nib = MsusY =500 GeV, T =200 GeV, (c) different squark masses 
niQ = 500 GeV,mt = 800 GeV,; nib = 200 GeV at zero temperature, and (d) different squark masses 
niQ = 500 GeV, mt = 800 GeV, mb = 200 GeV at T =200 GeV. 
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Figure 17: Effective temperature-dependent parameters Aj {i=l,...7) witli tlie one-loop thresliold 
and logaritlimic corrections at mz = 91.19 GeV, mi, = 3 GeV, rrit = 175 GeV, m\Y = 79.96 GeV, 
c/2 = 0.6517, gi = 0.3573, Gp = 1.174 • 10"^ GeV'^, MsusY = 500 GeV, mg = 500 GeV, nit = 800 
GeV, rub = 200 GeV, =2000 GeV, A = Xt + fi/tan(3, Xt = 700 GeV, tan/3 = 5, ht = I, = 0.1. 
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